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Abstract 

We consider three modifications of our involutive algorithm for computing Janet 
bases. These modifications are related to degree compatible monomial orders and 
specify selection strategies for non-multiplicative prolongations. By using the stan- 
dard data base of polynomial benchmarks for Grobner bases software we compare 
the modifications and confront them with Magma that implements Faugere's F4 
algorithm. 

1 Introduction 

In [1] we designed an algorithmic approach to computing Grobner bases based on the new 
notion of involutive monomial division. This notion generalizes the constructive ideas 
of French mathematician M.Janet [2] on partition of independent variables for partial 
differential equations (PDEs) into multiplicative and non-multiplicative. He applied this 
partition to complete systems of PDEs to involutivion. Janet's receipt of the separation 
of variables generates a particular involutive division which we called in [1] Janet division. 

Our algorithms [3, 4], which improve those in [1], shown rather good computational 
efficiency when specialized for Janet division. Moreover, for the present, none of other 
known involutive divisions demonstrated its computational superiority over Janet division. 
Though the improved algorithms allow a user to output a reduced Grobner basis without 
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extra reductions in the final involutive basis, it is computationally worthwhile to minimize 
Grobner redundancy in the intermediate and in final involutive basis. Recently [5], we 
succeeded in such an improvement of Janet division, and called the improved division as 
Janet-like. The last division, since it does not partition the variables into multiplicative 
and non-multiplicative, is not involutive. But even so, it is very close to Janet division and 
preserves all its computational merits. Quite often, e.g. for most of benchmarks in the 
standard data base [6, 7], there is no substantial computational superiority of Janet-like 
division over Janet divisions. But for other classes of polynomial problems such as toric 
ideals this superiority can be enormous [5]. 

Apart from improvement of the division, there is another important source of optimiza- 
tion in the involutive algorithms: selection of non-multiplicative prolongations. The last 
constructions, i.e., the products of intermediate polynomials and their non-multiplicative 
variables, in the case of their involutive head reducibility, play in the involutive algorithms 
the same role as S— polynomials play in Buchberger's algorithm [8]. In the latter case any 
S— polynomial can be selected at any step of the algorithm. This enormous arbitrariness, 
as well as that in the reduction strategy, is one of the main obstacles on the way of opti- 
mization of Buchberger's algorithm. Only after many years of research some heuristically 
good selection strategies for S— polynomials, such as "sugar" [9], were found. 

As to the involutive approach, the reduction sequence is uniquely defined and an 
admissible choice of a non-multiplicative prolongation is subject to certain restrictions [4]. 
Nevertheless, for examples large enough, one can choose from many possible prolongations. 
For example, in the 7th order cyclic root example [6, 7] at the intermediate algorithmic 
steps there arise several hundreds prolongations such that any of them can be chosen. 
By this reason it is important to investigate a heuristical efficiency of different selection 
strategies. 

In the given paper we present three different selection strategies which as we found are 
computationally good. In so going, we restrict ourselves with degree compatible monomial 
orders. In practice, this is a reasonable restriction. It is well-known that heuristically 
best way of computing a Grobner basis for an arbitrary order is to compute it, first, 
for a degree compatible order and then to convert the basis into the desirable one by the 
FGLM algorithm [10] or by the Grobner walk [11]. After some preliminary definitions and 
conventions (Sect. 2) we modify in Sect. 3 the involutive algorithm [4], as specialized for 
Janet division, in accordance to those selection strategies. In Sect. 4 we give experimental 
comparison of the three modifications of involutive algorithm on the benchmarks from [6, 
7]. Here we also show the corresponding timings for the implementation in Magma [12] 
of Faugere's F 4 algorithm [13]. We conclude in Sect. 5. 

2 Preliminaries 

Throughout this paper we shall use the following notations and conventions: 

N> is the set of nonnegative integers. 

X = {xi, . . . , x n } is the set of polynomial variables. 
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R = K[X] is a polynomial ring over a zero characteristic field K. 
Id(F) is the ideal in R generated by F C R. 

M = {x*i ■ ■ ■ x l ™ | ik G N>o, 1 < k < n} is the monoid of monomials in R. 

degj(w) is the degree of Xi in u G M. 

deg(w) = X/ILi degj(w) is the total degree of u. 

mindeg(t/) = min{ deg(wi), . . . , deg(ufc) } for U — {u\, . . . , u k } C M. 

>- is an admissible monomial order such that x x >- x 2 y ■ ■ ■ y x n and for u, v G M 

deg(w) > deg(w) =>- u y v. 
If monomial m divides monomial w and deg(w) < deg(f), i.e. u is a proper divisor of 

we shall write u \Z v. 
lm(/) and lt(/) are the leading monomial and the leading term of / G R \ {0}. 
lm(-F) is the leading monomial set for FcK \ {0}. 
J is Janet division. 

Mj(u, U) is the set of J— multiplicative variables for monomial u G U C M. 
NMj(u, U) is the set of J— non-multiplicative variables of monomial u G £7 C M. 
J(u, U) is the submonoid in M generated by the power products of variables in Mj(u, U). 
u G U (U C M is finite) is a Janet divisor of v G M if v = u ■ w, w G J(u, U). 
HNFj(f, F) is the J(anet)— head normal form of / G R modulo FcK. 
NFj(f, F) is the J(anet)— (full) normal form of / G R modulo F C R. 

The following definitions are taken from [1]. 

Definition 2.1. Janet division. Let {7 C M be a finite set. For each < i < n partition 
U into groups labeled by non-negative integers d , . . . , di (U = [0]). 

[d , d u . . . , di] ■= {u E U \ d = 0,d 1 = deg 1 (u), • • • , = deg^w)}. 

Indeterminate Xj is J (anet) -multiplicative for n G (7 if m 6 [do, • • • , di-i] and degj(-u) = 
max{deg i (v) | u G [d , • • • , di-i]}- 

Definition 2.2. [J —reduction). Given a monomial order a finite set F G R \ {0} of 
polynomials and a polynomial p G R \ {0}, we shall say that: 

(i) . p G R is J— reducible modulo f G F if p has a term t = au (a G K, u G M, a 7^ 0) 

such that lm(/) is J— divisor of w, that is, u = lm(/) ■ v where v G J(lm(/), \m(F)). 
It yields the J— reduction p ^ g := p — (a/lc(f)) f ■ v. 

(ii) . p is J— reducible modulo F if there is / G F such that p is J— reducible modulo /. 

(iii) . p is in the J— head normal form modulo F (p = HNFj(p, F)) if \m(p) has no 

J— divisors in lm(F). 

(iv) . p is in the J— normal form modulo F (p = NFj(p, F)) if p is not J— reducible 

modulo F . 
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Definition 2.3. (J— basis) A polynomial set F is called J — autoreduced if each term in 
every / G F has no J— divisors in lm(F) \ lm(/). A J— autoreduced set F is called a 
Janet basis of Id(F) if 

(V/ G F) (Vx G NMj(f, F)) [ NFj(f ■ x,F) — 0] . (1) 

A Janet basis G is called minimal if for any other Janet basis F of the same ideal the 
inclusion lm(G) C lm(F) holds. 

A Janet basis is a Grobner one [1], though generally not reduced. However, similarly 
to a reduced Grobner basis, a monic minimal Janet basis is uniquely defined by an ideal 
and a monomial order. In that follows we deal with minimal Janet bases only and omit 
the word "minimal". 

3 Modified Involutive Algorithm 

Our Janet division algorithm is given by (cf. [3, 4]): 



Algorithm: JanetBasis (F, -<) 



Input: F G K \ {0}, a finite set; -<, a degree compatible order 


Output: G, a Janet basis of ld(F) 


1 


choose / G F of the minimal deg(lm(/)) 


2 


G:={/} 


3 


Q:=F\G 


4 


do 


5 


h := 


6 


while Q ^ and /i = do 


7 


choose p E Q with minimal lm(p) w.r.t. y 


8 


Q:=Q\ M 


9 


h:=NFj(p,G) 


10 


od 


11 


if /i ^ then 


12 


for all {g G G lm(g) □ \m(h)} do 


13 


Q:=QUfo}; G:=G\{ fl } 


14 


od 


15 


G :=GU{h} 


16 


Q :=Qu{g-x\ g e G, x G iVMj(lm(p), lm(G)) } 


17 


fi 


18 


od while Q ^ 


19 


return G 
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In its improved form [4], at the initialization step, i.e., before starting the main loop 4-22 
and after its modification in the loop, the set Q is J— head reduced modulo G. Now, as 
the first modification of the above algorithm JanetBasis we shall use only partial head 
reduction of elements in Q as shown in the following algorithm: 



Algorithm: JanetBasis I (F, -<) 



Input: FsR \ {0}, a finite set; -<, a degree compatible order 


Output: G, a Janet basis of ld(F) 


1 


choose / G F of the minimal deg(lm(/)) 


2 


G:={f} 


3 


Q:=F\G 


4 


do 


5 


S := { q G Q \ deg(lm(g)) = mindeg(lm(Q)) } 


6 


P:=0; Q:=Q\S 


7 


for all s G S do 


8 


S:=S\{s}; p:=HNFj(s,G) 


9 


if p 7^ then 


10 


P := PU{p} 


11 


fi 


12 


od 


13 


while P ^ do 


14 


choose p E P with minimal lm(p) w.r.t. >- 


15 


P:=P\{p}l h:=NFj{p,G) 


16 


for all {# G G lm(g') □ lm(/i)} do 


17 


g:=guM; G:=G\{^} 


18 


od 


19 


G := GU{/i} 


20 


Q:=QU{g-x| 9 GG,xG iVMj(lm(p), lm(G)) } 


21 


od 


22 


od while Q ^ 


23 


return G 



In this algorithm at step 5 all the polynomials in Q of the minimal head degree are 
collected in set S and then are J— head reduced modulo G in the for loop 7-12 with the 
collection of nonzero head reduced polynomials in P at step 10. Then in the while loop 
13-21 the polynomial in P with the least leading term is sequentially selected at step 14 
of the loop and inserted in set G after its tail J— reduction of step 15. As well as in 
algorithm JanetBasis the displacement of some polynomials from G to Q done at step 
17 provides minimality of the output Janet basis [1]. 

Subalgorithms HNFj and NFj which are called at steps 8 and 15 compute the Janet 
head and the full normal forms, respectively, in accordance with Definition 2.2. 
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Apparently, the modifications done in algorithm JanetBasis I in comparison with algo- 
rithm JanetBasis do not violate its correctness. As well as for the latter algorithm, when 
the main do- while loop terminates, the polynomial set G satisfies the conditions (1) in 
Definition 2.3. 

We consider now another modification of algorithm JanetBasis shown in the form of 
algorithm JanetBasis II. 



Algorithm: JanetBasis II (F, -<) 



Input: F6l\ {0}, a finite set; -<, a degree compatible order 


Output: G, a Janet basis of ld(F) 


1 


choose / E F of the minimal deg(lm(/)) 


2 


G:={f} 


3 


Q:=F\G 


4 


do 


5 


S := { q E Q \ deg(lm(g)) = mindeg(lm(Q)) } 


6 


P:=0; Q:=Q\S 


7 


for all s E S do 


8 


S:=S\{s}; p:=NFj(s,G) 


9 


if p 7^ then 


10 


P := PU{p} 


11 


fi 


12 


od 


13 


P := Update(P, -<) 


14 


for all p E P do 


15 


for all {g E G \ lm(g) □ lm(p)} do 


16 


Q:=QU{g}; G:=G\{g} 


17 


od 


18 


G := G U {p} 


19 


Q:=QU{ 5 -x| 9 GG,xG NMj(\m(g),\m(G)) } 


20 


od 


21 


od while Q ^ 


22 


return G 



As distinct from algorithm JanetBasis I, the full J— normal form is computed at step 
8. Besides, at step 13 the polynomial set P whose elements are inserted in G in the for 
loop 14-20 is updated in accordance to the below subalgorithm Update. 

At steps 1 and 4 of the subalgorithm we indicate two different options for the choice 
of element p: with the highest or with the lowest leading monomial with respect to the 
order y. In our numerical experiments presented in the next section we apply these two 
different upgrade strategies when only one of the indicated choices (highest or lowest) is 
used in the whole run of the algorithm. 



6 



Algorithm: Update(P, >-) 



Input: FcK \ {0}, a finite set; >-, an order 


Output: H C R \ {0}, an updated input set 


1 


choose / £ P with the highest /lowest lm(/) w.r.t. >- 


2 


H :={/}; P:=P\{/} 


3 


while P ^ do 


4 


choose p G P with the highest /lowest lm(p) w.r.t. >- 


5 


P:=P\{p} 


6 


h:=NFj(p,H) 


7 


if /i ^ then 


8 


P :=ffU{/i} 


9 


fi 


10 


od 


11 


return P 



In subalgorithm Update an element / in the input polynomial set P (which is J— reduced 
modulo polynomial set G when the subalgorithm is invoked in the main algorithm) is 
chosen at the initialization step 1 with the highest or lowest leading term, depends on 
the selection strategy used. After that, in the first run of the while loop 3-10 the other 
polynomial p in P, if any, with the same leading monomial as that in / is J— reduced 
modulo /. In the case of nonzero reduction (when monic p is different from monic /) the 
normal form obtained is added to / at step 8 to be involved in the further reductions. 
Then, the processes of the selection and J— reduction of elements in P is repeated until 
P becomes empty. 

The above described modifications related to certain selection strategies for non- 
multiplicative prolongations are easily adapted to the improved version of involutive al- 
gorithm [4]. The improved version avoids useless repeated prolongations and applies the 
involutive analogues of Buchberger's criteria for detection of some superfluous reductions. 
Furthermore, it is straightforward to include the modifications into the Janet-like division 
algorithms [5]. 

4 Computer Experiments 

The improved version of algorithm JanetBasis I was implemented in C on as a part 
of package JB [3] whose version is also included in the library of Singular [14] and in 
C++ as a part of the open source software GINV [15]. The last software implements also 
algorithm JanetBasis II in its improved version and for both options in subalgorithm 
Update. For all that Ginv implements also Janet-like division [5]. 

We examined the three selection strategies of Sect. 3 by the standard data base of poly- 
nomial benchmarks [6, 7] and for degree-reverse-lexicographical monomial order. Some 
of the benchmarks are listed in the below table together with the timings they took for 
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computing Grobner bases. For comparison, we also included the timings of the last two 
versions of Magma [12] with the fastest Grobner bases module among all computer algebra 
systems. This is owing to implementation of the Faugere _F 4 algorithm which rests upon 
linear algebra for doing multiple reductions in contract to Buchnerger's algorithm [8] or 
our involutive algorithms which are relayed on the chains of elementary reductions. 

As we noticed in Introduction (Sect.l), the involutive algorithm in its improved form [4] 
can output a reduced Grobner basis as an internally fixed part of Janet basis without any 
extra reductions. By this reason all software included in the table output the same bases. 

It should be noted that GINV for can also use Janet-like division instead of Janet 
division whereas the package JB implements Janet division only. Generally, intermediate 
polynomial sets for Janet-like division are more compact then those for Janet division. 
There are whole classes of interesting multivariate polynomial problems, for example, 
binomial toric ideals [16] closely related to integer programming [17] for which Janet-like 
division leads to enormous gain in comparison with Janet division. However, as we already 
mentioned (Sect.l), the difference of two divisions does not manifest itself significantly 
for benchmarks in the table. 

The timings in the table were obtained on the following machines: 

JB: 2xOpteron-242 (1.6 Ghz) with 4Gb of RAM running under Gentoo Linux 2004.3 
with gcc-3.4.2 compiler. 

GINV: Turion-3400 (1.8 Ghz) with 2Gb of RAM running under Gentoo Linux 2005.1 
with gcc-3.4.4 compiler. 

Magma: dual processor Pentium III (1 Ghz) with 2 GB of RAM running under SuSE 
Linux 8.0 (kernel 2.4.18-64GB-SMP) with gcc-2.95.3 compiler. 

All timings in the table are given in seconds, and (*) shows that the example was not 
computed because of the memory overflow. 

The 2nd and 3th columns in the table show the results for algorithm JanetBasis 
I whereas the 4th and 5th columns shows those for algorithm JanetBasis II with the 
choice of the highest and lowest option in subalgorithm Update, respectively. 

One can see rather high stability of the involutive algorithm with respect to three 
variations used for the selection strategy. In addition, in all three cases the number of 
redistributions between G and Q was experimentally tested to be minimal. In so doing, we 
observed that the strategy in algorithm JanetBasis II with the lowest element choice in 
subalgorithm Update leads to a slightly more smooth growth of the intermediate memory 
needed than its counterpart with the highest element choice. The hcyclic8 example in the 
table explicitly demonstrates this observation. 

As to comparison with Magma, it clearly signals on superiority of the linear algebra 
based algorithm over our reduction strategy that uses the elementary reduction chains. 
In our future research plans there is also improvement of the involutive algorithm by doing 
reductions by means of linear algebra. 
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Benchmarking 



Example 


Algorithm I 


Algorithm I 


Algorithm II 


Algorithm II 


Magma 


Magma 




(JB) 


(GINV) 


high (GINV) 


low (GINV) 


V2.11-8 


V2. 12-17 


assur44 


10.35 


14.20 


6.33 


6.4 


4.56 


4.99 


butcher8 


1.06 


1.02 


0.38 


0.39 


4.68 


5.00 


chemequs 


0.67 


0.61 


0.57 


0.6 


12.80 


11.99 


chemkin 


17.83 


16.87 


10.95 


9.95 


32.34 


29.83 


cohn3 


76.72 


107.14 


30.21 


25.47 


37.73 


39.20 


cpdm5 


1.78 


1.57 


1.69 


1.68 


0.69 


0.70 


cyclic6 


0.12 


0.19 


0.14 


0.14 


0.09 


0.08 


cyclic7 


58.72 


60.94 


68.59 


65.28 


6.64 


7.08 


cyclic8 


12056.24 


14046.26 


5826.18 


4424.96 


235.73 


245.65 


dl 


8.77 


12.58 


1.99 


2.08 


28.49 


8.29 


desl8.3 


0.19 


0.18 


0.19 


0.19 


1.81 


1.89 


des22_24 


0.68 


0.62 


0.77 


0.79 


1.37 


1.46 


discret3 


23322.8 


20956.31 


12642.49 


13521.65 


33658.09 


19369.53 


dl 


270.17 


278.89 


80.77 


89.52 


14.57 


11.95 


eco8 


0.40 


0.44 


0.44 


0.46 


0.20 


0.20 


eco9 


3.22 


5.60 


4.99 


5.08 


1.25 


1.20 


ecolO 


52.56 


56.70 


65.71 


68.06 


7.07 


6.91 


ecoll 


765.98 


741.74 


718.53 


679.3 


62.33 


51.08 


extcyc5 


1.35 


1.53 


1.46 


1.37 


0.37 


0.38 


extcyc6 


324.70 


184.49 


276.06 


155.64 


45.36 


47.96 


extcyc7 










8242.00 


8492.13 


f744 


4.88 


7.71 


2.22 


2.68 


1.47 


1.38 


f855 


132.97 


139.79 


37.64 


38.45 


48.63 


37.06 


fabrice24 


108.52 


116.77 


8.2 


7.7 


9.45 


8.70 


filter9 


20.97 


5.76 


1.13 


1.6 


80.04 


56.67 


hairer2 


62.91 


108.17 


126.69 


125.43 


92.07 


85.86 


hairer3 


1.96 


0.92 


0.32 


1.4 


* 


* 


hcyclic7 


64.17 


53.87 


65.81 


73.0 


6.26 


6.76 


hcyclic8 


6024.97 


4316.59 


* 


7560.99 


229.70 


237.12 


hf744 


22.17 


8.58 


7.18 


11.26 


1.39 


1.32 


hf855 


2157.88 


534.08 


806.51 


988.38 


48.15 


36.69 


hietarintal 


0.77 


0.71 


0.38 


0.53 


2.63 


2.15 


il 


98.24 


122.36 


58.29 


58.21 


55.07 


42.35 


iliasl3 


1167.18 


5851.97 


3013.1 


2469.62 


336.21 


309.64 


ilias_k_2 


323.59 


669.68 


445.51 


270.21 


55.41 


54.71 


ilias_k_3 


452.32 


846.19 


1162.7 


622.14 


90.67 


89.97 


jcf26 


224.96 


211.24 


16.44 


14.65 


31.64 


25.59 


katsura7 


2.15 


1.77 


2.08 


1.98 


0.72 


0.79 


katsura8 


27.48 


24.66 


28.8 


27.09 


4.7 


5.06 


katsura9 


337.52 


294.59 


340.45 


311.98 


33.47 


34.87 


katsuralO 


4790.55 


4983.11 


7220.29 


6204.95 


287.38 


292.02 


kinl 


15.18 


20.32 


7.11 


7.11 


50.56 


45.33 


kotsireas 


6.33 


37.94 


4.93 


4.27 


3.45 


3.67 


1100116 


0.97 


1.29 


1.27 


1.29 


0.60 


0.62 


1100117 


28.87 


32.58 


37.52 


38.52 


4.93 


4.77 


1100118 


1552.26 


2292.84 


3322.62 


3152.57 


43.65 


42.80 


pinchonl 


10.37 


0.04 


0.01 


0.01 


4.09 


3.54 


rbpl 


210.94 


177.51 


173.8 


173.98 


38.33 


35.79 


rbpl24 


108.78 


116.78 


8.23 


7.7 


9.62 


8.74 


redcyc6 


0.16 


0.17 


0.13 


0.14 


0.10 


0.10 


redcyc7 


913.75 


1048.69 


48.19 


48.61 


5.73 


6.36 


redecolO 


18.51 


18.66 


23.91 


22.4 


2.33 


2.40 


redecoll 


178.32 


187.36 


253.34 


228.41 


14.56 


14.85 


redecol2 


1735.95 


2172.75 


4666.8 


3385.97 


101.51 


103.02 


reimer5 


0.22 


0.36 


0.34 


0.38 


0.74 


0.70 


reimer6 


9.69 


21.60 


24.19 


23.96 


42.13 


42.40 


reimer7 


719.37 


3808.91 


4756.4 


4314.12 


5216.53 


5032.73 


virasoro 


9.69 


8.90 


10.96 


10.68 


1.72 


1.77 
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5 Conclusion 



In this paper we experimentally investigated three different selection strategies for the in- 
volutive algorithm specialized for Janet division and observed its computational stability 
with respect to these variations in selection strategy. However, the problem of finding 
heuristically best selection strategies for the Janet division algorithms as well as for algo- 
rithms exploiting other involutive division is still open and is of practical importance. Our 
computer experimenting shows that the arbitrariness in selection of non-multiplicative 
prolongation at the intermediate steps of the algorithm is sharply grows with the number 
of variables and degree of the initial polynomials. 

By this reason the number of prolongations with the same head degree or even with 
the same leading monomial may achieve many hundreds and thousands for sufficiently 
large examples. That is why, searching for heuristically best strategies is so important for 
increasing computational efficiency of the involutive algorithmic methods. 
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